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Abstract 



' We prove a 2-terms Weyl formula for the counting function N{fi) of 

■ the spectrum of the Laplace operator in the Euclidean disk with a sharp 

remainder estimate O (/i^^'^). 

^ : Introduction 
cn 

^ ■ Let us denote by < Ai < A2 < ■ • ■ the eigenvalues for the Dirichlet Laplacian 

CN ■ of some bounded connected smooth domain X in the Euclidean plane. It has 

Tj- ! been shown by Ivrii [11] that the following 2-terms Weyl formula holds under 

some genericity assumption on the periodic orbits of the associated billiard ball 
problem: if Nx{fi) = #{j | A^- < /i^}, 



with R{n) = o{fi). Moreover, this result is quite optimal: Lazutkin and Terman 
[T3] showed that there is no 5 > so that an estimate R{fi) = O holds for 

all smooth convex domains. 

Our goal is to get an upper bound for R{fJ.) in the case of the Euclidean disk. 
Our main resu lt0 is: 

Theorem 1 

iVdisk(/i) = /iV4 -1^/2 + (^2/3) . 



*Institut Fourier, Unite mixte de recherche CNRS-UJF 5582, BP 74, 38402-Saint Martin 
d'Hcrcs Cedex (France); yves.colin-de-verdiere@ujf-grcnoble.fr 

^ After having completed this work, we learned from I. Polterovich that the same result has 
been announced in 1964 by N. V. Kuznecov and B. V. Fedosov in [H]. The method is similar 
to ours. We give here an independent complete derivation of the needed Van der Corput result 
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The proof is based on the exphcit expression of the eigenvalues as the squares 
of the zeros of the Bessel functions J„ as well as on some precise asymptotics of 
these zeros which goes back to Olver (see [15l |5]). This way, we have to study a 
lattice point problem in some domain with cusps. A rather general lattice point 
problem was studied by Van der Corput [16] , see also [101 El [T71 E] . In [1] , a similar 
method was used in order to get a good remainder estimate for some surfaces 
of revolution. Let us note also that the same remainder estimate holds for the 
integrable polygonal billiards like the rectangles or the equilateral triangles: this 
is a direct consequence of the explicit formula for the eigenvalues which reduces 
the question directly to a lattice point problem for which the Van der Corput's 
result applies. 



1 The spectrum of the unit disk 

We consider the spectrum of the Euclidean Laplacian Adisk = — <9| — in the 
unit disk in y with Dirichlet boundary conditions. As it is well known and can 
be checked by separation of variables, the eigenvalues of Adisk are the squares of 
the zeros of the Bessel functions J„, n eTj. Let us recall that 

U^) = ^ r e*(^^"^*-"*)rft , (1) 

and that we have the following identities Jn{—x) = (— l)"J„(x), J_„(— x) = 
J„(x). Let us denote by \n\ < Xi{n) < X2{n) < ■ ■ ■ < Xk{n) < ■ ■ ■ the positive 
zeros of J„. Then the spectrum of Adisk? with multiplicity, is given by 

o" = {xkin)'^ I n G Z, k = 1, ■ ■ ■ ,} . 



In order to describe the asymptotics of the zeros of Bessel functions, we in- 
troduce the domain D in defined by 

D = {{x,y) I — l<x<l, y < g{x), y > max(0, — x) } 

with ^ 

g{x) = — ^Vl — .x^ — X arccos . 

Let us define 71 = {{n,k — 1/4) | {n,k) G Z^} and S = {{x,y) \ y > 
max(0,— x)}. Let F : 5 — M be the function homogeneous of degree 1 which 
satisfies F = 1 on the graph of g. The spectrum of the disk is approximately 
given by {\j | j G N} ~ {F(m) | m G 7^ fl S}. More precisely, for n > 0, 
Xk{n) ~ -F(n, k — 1/4), while for n < 0, Xfc(n) ~ F(n, k + \n\ — 1/4). This will 
reduce our problem to a lattice point problem: 

iVdisk(/i) ~ NM = #{"^ G 7^ n /iD} . 
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y 




Figure 1: the domain D 



Remark 1 Let us note that D and IZ are invariant by the linear involution 
J{x,y) = {—x,y + x). This corresponds to the fact that Jn and J^n have the 
same zeros. 

In order to complete the argument, we will have to study the lattice point 
problem (Section [2]) and to show how close A^diskl/^) and iVD(/i) are (Section 
[3]). The first part uses the method of Hlawka, Herz and Randol [IHl [9l [17] for 
studying smooth lattice point problems and the second part is done using Olver's 
asymptotics for the zeros of Bessel functions: we re-derived it in [5] using the 
integral representation of the Bessel functions and the general theory of oscillatory 
integrals associated to versal unfoldings of singularities as explained in [H]. 

2 A lattice point problem with a cusp 

Let us denote by 7^ the lattice 7^ := {(n, A; - /3) | (n, A;) G Z^} with < /3 < 1. 
Let us consider a domain G C with a cusp: G = {{x,y) |0<x<l, < y < 
g (x)} with g (x) ~ a (1 — x)^^'^ with a > near x = 1. We consider the weighted 
lattice point problem defined by the counting function 

m={mi,m,2)(zfJ.GmZ 

with X ^ C*^(M) with x = 1 near and the support of x small enough. 
We have the following 2-term Weyl estimate: 

Theorem 2 Under the previous assumptions on G, (3 and x, we have 
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Corollary 1 If D is the domain defined in SectionU\ and /3 = 1/4, we have 

Proof of Corollary [I].- we decompose into 3 terms: one for each cusp and 
one inner term using an homogeneous partition of unity. The corollary follows 
from the previous Theorem for the parts near the cusps and from the classical 
estimates going back at least to Van der Corput [16] (see also [171 Ej) ioi the 
inner part. In order to use Van der Corput estimates O (/i^^^), we need to check 
the strict convexity, in fact the non vanishing of the curvature of the graph of g: 

this comes from the fact that g"{x) = {1 — x"^) ^ > 0. □ 
Proof of Theorem\^ let us denote by B{m,r) the Euclidean ball of center 
m and radius r. We can first replace xiu/^) by the smooth function Xoi^^v) = 
x{y/x){l — (f){x, y)) with (p G with support in the ball -8(0, min(/3, 1 — (3)) and 
= 1 near because there is no element of TZ in the support of 0. The smooth 
function xo is a classical symbol of degree 0: dldyXo{x,y) = 0{{l + \x\ + \y\)~^^~^^^). 

Let us give a positive function p G (M^) with Support (p) C {x^ + y"^ < 1} 
and /jg2 p{x, y)dxdy = 1, define Pe = p {-/s) / with e = p~^^^, and consider 

^.^(/^)= E {xolGi.^Pe)M (2) 

men 

where 

Gte = {i^^y) \ 0<x<p, 0<y<pgix/p) + 2e} , 




For each m ^ p.G with m eTZ, B (m, e) fl G'^ = 0, hence {1q- Pe){'m) = 
while \/{x,y) G M?, < (1^-^ * Pe){x^y) ^ 1- Similarly, for each m G pG fl 71, 
B {m,e) C G^j, and (1^+^ *Pe)("^) = 1- Hence, 

ip) <Nip)< iV+ (p) . 
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We will apply Poisson summation formula and use estimates on the Fourier trans- 
form of . Let us denote by 

the Fourier transforms of Xo-times the characteristic function of G^^. 
The Poisson summation formula applied to the sum ([2]) gives 



/ Xo{x,y)l(.± {x,y)dxdy+ V p (27r£ (p, g)) $J (27rp, 27rg) e" 



N^{fi)= / xo{x,y)lG^ {x,y)dxdy+ ^ n (o^^ (r. <h± (o^^ o^n\ ^-■^^^Pi 

(p,g)ez2\o 

(3) 

We need to evaluate We use Green- Riemann formula in order to eret 

integrals on the boundaries. We have the following formulas: 

Lemma 1 // xo is a smooth classical symbol of degree and 



Xo{x,y) + -dyXo{x,y) - \dyyX^{x,y) \ ^^''^^y'^Mx , 



I 

7] \ 1] 

then 

da = Xoix, y)e''^^^~^y'^^dx A dy + 0(?7~^)xi(a;, y)dx A dy , 

where Xi{x,y) G L^{dxdy). 
A similar results holds for 

^ = ^ (^Xo{x, y) + jd,^Xo{.x, y) - ^d^^xo{x, y)^ e^^^^+^^^^dy 

If |?7| < C\^\, we use 

e^i'^^+y'^^Xoix, y)dx Ady = d/3 + 0{l/^^)x2dx A dy , 
while, if |,^| < C|?7|, we use 

eiK+?/'?)^Q(2;, y)dx A dy = da + 0{l/ri^)xidx A dy . 
We have to estimate the integrals 

e^ii-+^y)xo {x,y)dy, 



where v = rj/^ is bounded (and similar integrals with xo replaced by the deriva- 
tives of Xo which are symbols of < degrees) and 

act. 



where u = ^/t] is bounded. We use the upper bounds given in Appendix A for 
the different parts of the boundaries, using the parametrization of the graph of 
y = g{x) by x{t) = 1 — t^f{t), y(t) = for < t < to- For example, the main 
part of the integral on the curved part of dG^^ of a is 



I 

V Jo 



to which we apply estimate given in Lemma [HI This gives: 
Lemma 2 The following estimates hold: 
m 

(0, 0) = /i^ / X (y/x) dxdy + O {fi'/') 

JG 

• For 1 < IpI < C\q\, 

$J,(27rp,27rg) = 0' 



For p = 0, q ^ 0, 



• For \q\ < C\p\, 

^te (2^P' 27rg) = O ( J" + ) 



Let us prove for example the estimate of (0, 27rg). The corresponding integral 
on the boundary splits into 2 parts ..dx — jj" ..dt. The first part gives the first 
term. The second part gives, up to constants, J = q~^ exp{27iiqnt^)xo{lJ'X{t) , ny{t)) j^x' {t)dt 
which is bounded by 0(g~^/i(g/i)~^/^) using the estimates of Lemma El 
We need also the classical formula: 

Lemma 3 For < /3 < 1, we have 

Theorem [2] follows then from the previous Lemmas and simple evaluations of 
the sums in the Poisson summation formula ([3]); using the fact that the Fourier 
transform of p is rapidly decaying, we need the bounds: 
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Lemma 4 We have: 

(p,(?)ez2\o 

J2 b?r^(i+/i-^/'ii(p,^)ii)-^ = o((iog/xr) , 

i<H<c^kl 

/^-'/' E br^(i+/i-^/'ii(P,g)ii)-"' = o(i) . 

i<kl<cbl 

Let us check the first upper bound, the others are similar. The first sum is 
bounded by 

Cf^^ E \\{p,q)\\-'^'{l + f^~'^'\\{p,q)\\r'^' 
which is of the same order as the integral 



-oo 

1 / 



rdr 



r3/2(i + /i-i/3r)^ ■ 

□ 



3 Spectrum of the disk as a lattice point prob- 
lem 

Our goal is to prove the following result: 
Theorem 3 

iVdi.k(/i) =iVD(/^) + 0(/i2/3) . 

This will complete the proof of Theorem [H 
Proof. - 

The estimate splits into 3 parts: the inner part and the 2 boundary 
parts. We choose a function x ^ — 1,1[, [0,1]) which is = 1 

in some large interval [— 1 + c, 1 — c] and split the two numbers to 
compare as 

iVdisk(/i) = iVd\sk(/^) + nLM + ivik(/i) = 

= E X{k/n)+ E i^-x{k/n))+ E i^-xik/n)), 

Xk{'n.)<li i^yO, Xk{n)<ii n<0, xi;{n)<fi 



and 

NM = N'M + N'M + N'M 



X{k/n)+ Yl {l-x{k/n))+ J] {l-x{k/n)) 



E 

(n,fc+max{0,-n)-l/4)G/iZ) n>0, (■n,k-l/4)£fiD n<0, {n,k+\n\-l/4)efiD 



We will compare the first terms (the inner parts) in both decomposi- 
tions and the second terms (the boundary parts). The third ones are 
similar to the second ones. 

The inner part: The zeros Xk{n) are given uniformly in any domain 
Xkiji) > (1 + c)\n\ with c > 0, by 



Xk[n] 



F{n, k-\) + 0(1/(1 + A; + n)) if n > 

F{n, k+\n\-\) + 0(1/(1 + k+ \n\)) if n < 



This is a consequence of the stationary phase expansion applied to 
the integral representations ([1]) of Bessel functions (see [5]). Using 
the fact that when F{n,k + max(0, —n) — 1/4) is close to fi, \n\ + k 
is of the same order as we get 

K {^'-^)^ ^disk(/^) <Nh(^fi + ^ 

It follows then from the Van der Corput's remainder estimate O [fi'^^^) 
for the smooth strictly convex lattice point problems that 

nLM - N'M = o i^.'^') . 

The boundary parts: due to the fact that the zeros of J„ and J_„ 
are the same, we discuss only the case n > 0. We are in a domain 
where Xk{n) < (l + 0)n. Let us denote by tk the k-th zero of the Airy 
function, we have (see pQ) 



tk 



2^(*-i 

2 V 4 



2/3 

+ eik) , 



with e{k) = O [k ^). From 0, we have the following equation for the 
zeros X of the Bessel function J„(x) : 

Ai (n)) + b (n, x) \x |-^/3Ai'(|x|2/V(n)) = , 

where m = (n/x) — 1 < 0, p is a smooth germ of odd diffeomorphism 
of (M, 0) (with p' > 0) and b {u, x) is a smooth symbol of degree in x. 
We deduce, using the implicit function theorem and the asymptotics 
of the Airy function, that 

Xk{n) = n{l + ip { ) ) + ri{n, k) 



n 



with ri{n,k) = 0(n~^), ip smooth, ip{0) = 0, ip'{0) > 0. This asymp- 
totics is due to Olver |T5]. If (1 + C)n > k > {1 + c)n > with 
< c < C , this asymptotics matches with the inner asymptotics via 
the asymptotics of the tkS for large /c's. 
Let 

Nk (/i) := k - 1/4) efiD\l<k<Cn} 

and 

K := <fi\l<k<Cn} . 

We have the 
Lemma 5 

liVfc (/i) - iV^ I < AT, (^^ + - AT, (/^ - ^) + Cfi'/'k-'/-' . 

By summing the estimate of the previous Lemma w.r. to k and using 
the 2-terms asymptotics of Nl){fi), we get 

nLM - N'M = o if,'/') . 

□ 

Proof of Lemma\^- Let us write F(x,y) := x{l + il)i{y'^/^ /x^^'^))] we have 
Arfc(/i) = A; - 1/4) < /i | 1 < A; < Cn} and (/i) = #{F(n, A; - 1/4 + 

e{k)) < fx + ri{k,n) | 1 < A; < Cn}, with e{k) = 0{l/k) and r/(A;,n) = 0(l/n). 
We have, in the range < y < cx with c small enough, < a < S^^F < h and 
= 0(x^/^i/~^/^). The Lemma follows by estimating the cardinal of the sets 

Ak:={n\ ii< F{n, k - 1/4) < /x + r/(A;, n) + Ce{k)n^l^k-^l'^} 

and 

A'^ := {ri I /i - n) - CE{k)n^l'^k~^l^ < F{n, A; - 1/4) < /i} . 

We use the fact that if F{n, — 1/4) is close to fi then n fi. We have Ak C 
5fc U Ck with 5fc := {n I /i < F(n, A; - 1/4) < /i + 0(l//x)} and := {n | /i + 
?7(A;, n) < F(n, A: — 1/4) </i + ?7(A:, n) + Ce(A;)n-'^/'^A;~-'^/^}. Using the estimate on 
d^F, we have = O (ni/^A:-^/^) 

4 Conclusion and problems 

It would be nice to get similar estimates for other integrable billiards like a 
circular annulus. The case of ellipse is more difficult and is due to Emile Mathieu 
[11]: the problem is with the unstable periodic geodesic (the larger diameter). 
We know now a good approximation of the associated eigenvalues thanks to my 
works with Bernard Parisse and San Vii Ngoc ([6l [7]). 
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Appendix A: Estimation of some integrals 



We need to get estimates of various integrals corresponding to part of the bound- 
ary of the domains D"!^^ to be defined in Section O 

Let us first recall the following stationary phase estimate: 

Lemma 6 Let f E C°°([a,6],C) and G C°°([a, 6], M) so that has only non 
degenerate critical points, then, if 




we have I{r) = O yr ^j. If (p depends smoothly on some parameter fi so that 
the non degeneracy assumption holds for fi = hq, the same conclusion is true 
uniformly in some interval |/i — /io| < c with c small enough. 

The curved part 

These integrals come when evaluating integrals on the curved part of the domains 

Lemma 7 Let us consider the integral 

hAr)= r e-(''''-''-f^')k'gm 
Jo 

with g e C^(]R), / E C°°(M,M) with /(O) ^ and \u\ < Uq < oo with z/q small 
enough, then, as r — )■ oo, /c,a(t) = 0(r~2) if \uq\ is small enough. 

This is easy using the stationary phase Lemma El 

Lemma 8 Let us consider the integral 

IcAr) = / e-i''-''''f^')Hg{t)dt 
Jo 

with g G C^(M), / G C°°(]R, M), and |z/| < z/q < oo with vq small enough, then, 
as T ^ oo, IcA^) = 0(r~2). 

This is more difficult because the critical point t = is degenerate for u = 0. We 
need a 

Definition 1 A smooth function f{t, a) is in S'' if all t-derivatives are hounded 
near t = oo by 0(t^) uniformly in a. 

Proof. - 
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We will first prove the Lemma for / = 1. Let us put jj, = ut^^^. We 
consider 2 cases: 

• < 1: let us make the change t = wt~^^^, we get 

Jo 

The critical points of the phase are and 2/i/3. We split the inte- 
gral into 2 parts with a smooth partition of unity 1 = h + {l — h) 
with h G C^([0, Wo[) and h = 1 on [0, 1]. The part containing the 
critical points is O (1) using an uniform bound for the integrand. 
For the other part, we introduce L = 1/ (3w^ — 2/xw) d/dw and 
integrate by part several times using the formal transpose *L of L 
and the fact that wg [wt~^^^) E S^, so that {*L)^ [wg {wt~^/^)) G 
^i-2Af_ If ]V > 2, this gives a function which is in L}{\wq, +00 [) 
uniformly in yU,. 

• l/^l > 1: we put t = ua and get 

roo 

I^j,{r) = / e*^'("'-'^')a(? [ya) da . 
Jo 

We split the integral smoothly and get for the part containing 
the critical points O (u'^ / fi^^'^^ = O (r~^/^). For the other part 
we use K = {l/3w^ - 2w) d/dw, : S^~^ and ag {ua) G 

S^. We pick a factor for each integration by parts and get 

It is clear enough that the proof still works if / is not constant. 



The linear parts 

These integrals come when evaluating integrals on the linear parts of the domains 



Lemma 9 We have, for \r]\ = 0{\^\), h {^,v) = f Jo' e'^^dy = O {e/\^\). 
For^ = 0{\ri\), zf 

V Jo 

then, for ^ ^ 0, 4 (e, v) = O {l/\^v\) M ^ ^ and I, (0, r/) = ^. 
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